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Abstract 

in 

We present a phase-space analysis of cosmology containing multiple scalar fields with a 
positive or negative cross- coupling exponential potential. We show that there exist power- 
law kinetic-potential-scaling solutions for a sufficiently flat positive potential or for a steep 
■^j- ' negative potential. The former is the unique late-time attractor, but it is difficult to yield 
assisted inflation. The later is never stable in an expanding universe. Moreover, for a steep 
negative potential there exists a kinetic-dominated regime in which each solution is a late- 
time attractor. In the presence of ordinary matter these scaling solutions with a negative 
cross-coupling potential are found unstable. We briefly discuss the physical consequences 
of these results. 
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1 Introduction 



Scalar field cosmological models are of great importance in modern cosmology. The dark 
energy is attributed to the dynamics of a scalar field, which convincingly realizes the goal 
of explaining current accelerating expansion of universe generically using only attractor so- 
lutions [L . Furthermore a scalar field can drive an accelerated expansion and thus provides 
possible models for cosmological inflation in the early universe j2j. In particular, there have 
been a number of studies of spatially homogeneous scalar field cosmological models with an 
exponential potential. They are already known to have interesting properties; for example, 
if one has a universe containing a perfect fluid and such a scalar field, then for a wide range 
of parameters the scalar field mimics the perfect fluid, adopting its equation of state [3]. 
These scaling solutions are attractors at late times [I]. The inflation models and other 
cosmological consequences of multiple scalar fields have also been considered jSJE]- 

The scale-invariant form makes the exponential potential particularly simple to study 
analytically. There are well-known exact solutions corresponding to power-law solutions 
for the cosmological scale factor a oc t p in a spatially flat Friedmann-Robertson- Walker 
(FRW) model [7j. More generally the coupled Einstein-Klein-Gordon equations for a single 
field can be reduced to a one-dimensional system which makes it particularly suitable for 
a qualitative analysis [H] E|. Recently, adopting a system of dimensionless dynamical vari- 
ables (TOj, the cosmological scaling solutions with positive or negative exponentials have 
been studied JIT] . In general there are many scalar fields with exponential potentials in su- 
pergravity, superstring and the generalized Einstein theories, thus multiple potentials may 
be more interesting. In the previous paper [T2j, We studied the stability of cosmological 
scaling solutions in an expanding universe model with multiple scalar fields with positive or 
negative exponential potentials. A phase-space analysis of the spatially flat FRW models 
shows that there exist cosmological scaling solutions which are the unique late-time attrac- 
tors and successful inflationary solutions driven by multiple scalar fields with a wide range 
of each potential slope parameter A. It is assumed that there is no direct coupling between 
potentials. Multiple cross-coupling exponential potentials arise in many occasions, for in- 
stance, from compactifications of vacuum Einstein gravity on product spaces [IB]. Indeed 
they are a natural outcome of the compactification of higher dimensional theories down to 
3+1 dimensions. With this in mind it is worth investigating such potential in a bit more 
detail. 

In this paper, we first study a system of dimensionless dynamical variables of multiple 
scalar fields with a positive or negative cross-coupling exponential potential. We obtain the 
scaling solutions and analyze their stability. There still exist cosmological scaling solutions 
which are the unique late-time attractors. In this model we then introduce a barotropic 
fluid to the system. We discuss the physical consequences of these results. 
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2 Cross-coupling Exponential Potential 



We consider n scalar fields (pi with a cross-coupling potential 

V = Vq exp I - £ > 



1=1 



where k 2 = 8ttGn is the gravitational coupling and A; are dimensionless constants charac- 
terising the slope of the potential. Further we assume all A, > since we can make them 
positive through 0; — > — <pi if some of them are negative. The evolution equation of each 
scalar field for a spatially flat FRW model with Hubble parameter H is 

4>i + SHfa - X iK V = 0, (2) 

subject to the Friedmann constraint 

H 2 = £ (f Ul + V] . (3) 



3 V^2 



Defining (n + 1) dimensionless variables 



Xi = -7=77, y = J J , (4) 



the evolution equations can be written as an autonomous system: 

^3 



n 

3 Xl (i-E^I ± A n/y> ( 5 ) 



/3 



where a prime denotes a derivative with respect to the logarithm of the scalar factor, 
N = In a, and the constraint equation (J3j) becomes 

E^ 2 ±y 2 = i- (7) 

i=i 

Throughout this paper we will use upper /lower signs to denote the two distinct cases of 
±Vo > 0. xf measures the contribution to the expansion due to the field's kinetic energy 
density, while ±y 2 represents the contribution of the potential energy. We will restrict 
our discussion of the existence and stability of critical points to expanding universes with 
H > 0, i.e., y > 0. Critical points correspond to fixed points where x\ = and y' = 0, and 
there are self-similar solutions with 

tt n 

4 = -3 £4 (8) 
n i=l 
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This corresponds to an expanding universe with a scale factor a(t) given by a oc t p , where 



P = o V n T 2 • ( 9 ) 

The system (jSJ) and © has at most one n-dimensional sphere S embedded in (n + 1)- 
dimensional phase-space corresponding to kinetic-dominated solutions, and a fixed point 
A, which is a kinetic-potential-scaling solution listed in Table 1. 

In order to study the stability of the critical points, using the Friedmann constraint 
equation (J7j) we first reduce Eqs.(jHJ and (JBJ) to n independent equations 



x\ = [XiJ-SxA i-$> 2 ■ (io) 




Substituting linear perturbations X{ — > Xj + 5xj about the critical points into Eqs. flTUj) , to 
first-order in the perturbations, gives equations of motion 

bx\ = -2 - 3xi j £( x i fa i) - 3 f 1 - £ x ? j ^ > (11) 

which yield n eigenvalues m.j. Stability requires the real part of all eigenvalues being 
negative. 

S: J27=i x l = 1, 2/ = 0. These kinetic-dominated solutions always exist for any form 
of the potential, which are equivalent to stiff-fluid dominated evolution with a oc 

£ l/3 

irrespective of the nature of the potential. Then Eqs.()ll|) become 




which yield n eigenvalues: one of them, say mi, does not vanish, mi = — y/6{Y%=i AjXj — a/6); 
the remains of them vanish. Thus the solutions are marginally stable for Y7i=i{^i x i) > 
a/6. For the special case Aj = A, using the constraint equation (JJJ) we find a/6/ (nX) < 
Y^=i x i/ n — Ol£=i x \ I 71 ) 1 — l/y/n. That is, if each scalar field has an identical-slope 
potential, there exist stable points only for A 2 > 6/n. 



A: X{ = y = y±(l — | Zir=i A 2 )- The potential-kinetic-scaling solution exists for 
sufficiently flat Y%=i A 2 < 6 positive potentials or steep Yh=i A 2 > 6 negative potentials. 
The power-law exponent, p = ^> 2 x 2 , depends on parameter Aj. From Eqs. ljllj) we find the 
eigenvalues 



rrii 
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Label 




y 


Existence 


Stability 


S 







all Aj 


stable Y,?=i(Kxi) > 


A+ 


A, 
V6 




E?=i A 2 <6, V>0 


stable 


A_ 


Aj 
V6 


\/(E?=if-i) 


E? =1 A? > 6, V < 


unstable 



Table 1: The properties of the critical points in a spatially flat FRW universe containing n 
scalar fields with the cross-coupling exponential potential. 



Thus the scaling solution is always stable when this point exists for a positive potential, 
which corresponds to the power-law inflation in an expanding universe when Y^=\ A? < 2. 
However, this solution is unstable for a negative potential. 

The different regions of Aj parameter space lead to different qualitative evolution. As 
an example we consider the cosmologies containing n scalar fields with the cross-coupling 
potential Aj = A. For the sufficiently flat (A 2 < 6/n) positive potential, these kinetic- 
dominated solutions are unstable and the kinetic-potential-scaling solution is the stable late- 
time attractor. Hence generic solutions start in the former and approach the later at late 
times. For the steep (A 2 > 6/n) positive potential, there exists a stable kinetic-dominated 
regime, in which each points are the late-time attractors. Hence generic solutions start in 
kinetic-dominated solution and approach the stable regime. For the flat sufficiently (A 2 < 
6/n) negative potential, only these kinetic-dominated solutions exist which are unstable 
scaling solutions. For the steep (A 2 > 6/n) negative potential, the kinetic-potential-scaling 
solution is unstable and there exists a stable kinetic-dominated regime. Hence generic 
solutions start in a kinetic-dominated regime or the kinetic-potential-scaling solution and 
approach the stable kinetic-dominated regime at late times. 



3 Plus a Barotropic Fluid 

We now consider multiple scalar fields with the cross-coupling potential (0) evolving in 
a spatially flat FRW universe containing a fluid with barotropic equation of state P 1 = 
(7 — l)p 7 , where 7 is a constant, < 7 < 2, such as radiation (7 = 4/3) or dust (7 = 1). 
The evolution equation for the barotropic fluid is 

p 7 = -3#(p 7 + P 7 ), (12) 

subject to the Fridemann constraint 

H 2 = Y (t^i+V + P,)- (13) 
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We define another dimensionless variable z = The evolution equations (J2J) and (|T2*)| 

can then be written as an autonomous system: 

* = -Sx^l-E^-pjiA,^ 2 , (14) 

n 



v = y\zj:^ + ^ 2 -Jl±^i), (is) 

i=l ' i=l ) 

- 7 + 2X>- +7^1, (16) 



3 

—z 
2 



i=i 

and the constraint equation becomes 



$>? ± j, 2 + = 1. (17) 



i=l 

Critical points correspond to fixed points where x\ — 0, y' = and z' = 0, and there are 
self-similar solutions with 

4 = -3E^-y- 2 - (18) 

i=l 

This corresponds to an expanding universe with a scale factor a(t) given by a oc t p , where 

2 

6 E i= i atf + 3iz 2 

The system ()14j ) -([16 |) has at most one n-dimensional sphere S embedded in (n + 2)- 
dimensional phase-space corresponding to kinetic-dominated solutions, a fixed point A 
which is a kinetic-potential-scaling solution, a fixed point B which is a fluid-dominated 
solution, and a fixed point C which is a fluid-potential-kinetic-scaling solution listed in 
Table 2. 

! x 2 — 1, y — 0, z — 0. These kinetic-dominated solutions always exist for any 
form of the potential, which are equivalent to stiff-fluid dominated evolution with a oc t 1 / 3 
irrespective of the nature of the potential. The linearization of system (|14j) - (jl6j) about these 
fixed points yields 



bx\ = -2 ( A^/~ - 3xi J ^(xjSxj), 



5z' = |(2- 7 )fo, 

which indicate that the solutions are marginally stable for Y^i=i{Kxi) > and a stiff fluid 
7 = 2. 
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A: Xi 



±(l - | E"=i a 4 2 . 



0. The potential-kinetic-scaling solution 



exists for sufficiently flat J27=i M < 6 positive potentials or steep E£=i Af > 6 negative 



potentials. The power-law exponent, p 



depends on the slope of the potential. 



V" , A 2 • 

The linearization of system (|14|) - (jl6|) about this critical point yields (n + 1) eigenvalues 



rrii 




Thus the scaling solution is stable for a positive potential with Yh=i Af < 37, which corre- 
sponds to the power-law inflation in an expanding universe when ElLi A? < 2. 

B: Xi = 0, y = 0, z = 1. The fluid-dominated solution exists for any form of the 
potential, corresponding to a power-law solution with p = 2/3^. 

5x[ = —35xi + (37 — y/6\i)8z, 
5z = 3^5z, 



which indicate that the solution is never stable. 
; A, 



C '. Xi 



scaling solution exists for a positive potential with E£=i A 2 > 37- The power-law exponent, 
p = 2/37, is identical to that of the fluid-dominated solution, depends only on the barotropic 
index 7 and is independent of the slope A» of the potential. The linearization of system 
(fTl j) -([TE |) about the fixed point yields 



3 (2-7)7 

2 > ™ ' A 2 ' 



Z 



3 7 



The fluid-potential-kinetic- 



Sx' 



by' 



3(2 - i)xt Yfritei) ~ 3 ( (1 - ?)(1 - E *i) + |^ 



+(\/6Ai 



H 2 -i)yJ2( x 3 Sx 3 



+ ?(2 -7)E^ 



j'=i 



v ^ i= i 

/3^ n , 3 7 9 7 2 



5y, 



which yield (n+l) eigenvalues 



mi 



m 2 



3(2-7) 
4 

3(2-7) 
4 



87(1X1 A| - 3 7 ) 
\T E? =1 A 2 (2- 7 ) 



87(E7=i Af - 3 7 ) 

\r Er =1 Af(2- 7 ) 



6 





X% 


y 


z 


Existence 


stability 


S 


TH A = i 








all Aj,7 


stable 7 = 2, 
E"(AiXj) > y/6 


A + 


A, 

Ve 


\/(i-E?$) 







stable (E™ Af < 3 7 ) 


A_ 


A, 


\/(E?$-i) 





e?a?>6 


unstable (V < 0) 


B 








1 


all Aj,7 


unstable 


C 


[3 7Ai 

V 2 y n a 2 


Iz (2-7)7 


V 1 A 


E?A? >3 7 


stable (V > 0) 



Table 2: The properties of the critical points in a spatially flat FRW universe containing n 
scalar fields with the cross-coupling exponential potential plus a barotropic fluid. 



m 3 = . . . = mn+1 = (2 - 7). 

Thus the scaling solution is stable for a positive potential with E™=i Af > 37. 

The different regions in the (7, Aj) parameter space lead to different qualitative evolu- 
tion. For the sufficiently flat (X"=i Af < 37) positive potential, S, A and B exist. Point 
A is the stable late-time attractor. Hence generic solutions begin in a kinetic-dominated 
regime or at the fluid-dominated solution and approach the kinetic-potential-scaling so- 
lution at late times. For the intermediate (37 < J27=i A? < 6) positive potential, all 
critical points exist. Point C is the stable late-time attractor. Hence generic solutions 
start in a kinetic-dominated regime, at the kinetic-potential-scaling solution or at the 
fluid-dominated solution and approach the stable fluid-kinetic-potential-scaling solution. 
For the steep (E"=i Af > 6) positive potential, S, B and C exist. Point C is the sta- 
ble late-time attractor. Hence generic solutions start in a kinetic-dominated regime or at 
the fluid-dominated solution and approach the stable fluid-kinetic-potential-scaling solu- 
tion. For the sufficiently flat (EILi Af < 37) negative potential, the kinetic-dominated 
solution S and the fluid-dominated solution B exist, which are unstable. For the inter- 
mediate (37 < Yh=i Af < 6) negative potential, the kinetic-dominated solution S and the 
fluid-dominated solution B exist, which are unstable. For the steep (J27=i Af > 6) negative 
potential, S, A and B exist. Point A is the stable late-time attractor. Hence generic solu- 
tions start in a kinetic-dominated regime or at the fluid-dominated solution and approach 
the stable kinetic-potential-scaling solution at late times. 
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4 Conclusions and Discussions 



We have presented a phase-space analysis of the evolution for a spatially flat FRW universe 
containing n scalar fields with a positive or negative cross-coupling exponential potential. 
In particular, for the Aj = A case, we find that in the expanding universe model with a 
sufficiently flat (A 2 < 6/n) positive cross-coupling potential the only power-law kinetic- 
potential-scaling solution is the late-time attractor. It is more difficult to obtain assisted 
inflation in such models since the fields with cross-coupling exponential potential tend to 
conspire to act against one another rather than assist each other. However, steep (A 2 > 6/n) 
negative cross-coupling potential has kinetic-dominated solutions with a oc some of 
which are the late-time attractors. It can be known that the kinetic energy of each field 
tends to be equal via their effect on the expansion at late times. 

Then we have extended the phase-space analysis of the evolution to a realistic universe 
model with a barotropic fluid plus n scalar fields with a positive or negative cross-coupling 
exponential potential. We have shown that for the sufficiently flat (X^=i ^? < ^7) positive 
cross-coupling potential, the kinetic-potential-scaling solution is the stable late-time attrac- 
tor. The energy density of the scalar fields dominates at late times. Moreover, for the steep 
(J27=i ^1 > 6) positive cross-coupling potential, the fluid-kinetic-potential-scaling solution 
is the stable late-time attractor. However, a negative cross-coupling potential has no stable 
scaling solutions. 
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